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PRESSURE AND FLOW- FIELD STUDY AT MACH NUMBER 8 


OF FLOW SEPARATION ON A FLAT PLATE WITH 
DEFLECTED TRAILING- EDGE FLAP 

By Charles B. Johnson 
Langley Research Center 

SUMMARY 

An experimental investigation was made of the flow separation on a flat plate model 
with a short trailing-edge flap deflected at angles of 10°, 20°, and 30° relative to the plate 
surface. These tests were conducted at a nominal free-stream Mach number of 8 and the 
nominal free-stream unit Reynolds number per foot was varied from 0.22 x 10® to 
10.9 x 10® (per meter from 0.72 x 10® to 35.8 x 106). Pressure measurements and 
schlieren studies were made for wall-to-total temperature ratios of 0.14, 0.43, and 0.74. 
Surface oil-flow studies were made at a wall-to-total temperature ratio of 0.43. Local 
similarity boundary -layer calculations were made upstream of the interaction region. 
Properties of the interaction and separated flow region were calculated with the boundary- 
layer— shock-wave interaction theory of Lees and Reeves (AIAA Journal, Nov. 1964). The 
upstream flat-plate boundary-layer theory was joined to the Lees and Reeves theory at 
the beginning of the interaction region by matching the boundary -layer momentum 
thicknesses. 

The results showed good agreement between the experimental and calculated pres- 
sures for wall-to-total temperature ratios of 0.74 (adiabatic) and 0.43 (room temperature). 
The oil-flow study showed the variation in the extent of separation with a change in unit 
Reynolds number for both laminar separation and for transitional separation. Wall 
cooling, for transitional separation, reduced the extent of separation. The maximum 
pressures measured on any particular flap varied considerably with a change in unit 
Reynolds number. The maximum pressures measured on the 30° flap at a unit Reynolds 
number per foot of 2.65 x 10® (per meter of 8.69 x 10®) indicated a pressure rise of 
approximately 58 percent over the oblique shock reattachment value. 

INTRODUCTION 

Flow separation is a common fluid mechanical phenomenon occurring on many con- 
figurations over a range of conditions from subsonic to hypersonic flow. Many investiga- 
tions, both theoretical and experimental, have been made of varied configurations and 
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conditions causing flow separation. Recent reviews of flow- separation research are 
given in references 1 to 3. As the free-stream Mach number increases into the hyper- 
sonic flow regime, laminar flow generally becomes more prevalent; therefore, increased 
attention has been focused on laminar separation and the associated laminar plateau pres- 
sure rise. Typical examples of experimental investigations of flow separation on flat 
plates and axisymmetric bodies are found in references 4 to 23. 

The purpose of the present paper is to present a detailed experimental investigation 
and a comparison with theory of the effect of unit Reynolds number, flap angle, and wall- 
to-total temperature ratio on the surface static pressure and the flow field for a flat plate 
model with a short trailing- edge flap. The tests herein were conducted for wall-to-total 
temperature ratios of 0.14, 0.43, and 0.74, flap angles of 10°, 20°, and 30°, and unit 
Reynolds numbers per foot ranging from 0.22 X 10® to 10.90 x 10® (per meter from 
0.72 x 10® to 35.76 x 10®). Pressure measurements and schlieren studies were made for 
all three wall-to-total temperature ratios. Previous tests conducted with a similar con- 
figuration have been made over a narrower range of unit Reynolds numbers and have not 
considered the effect of wall temperature. (See refs. 6, 7, 15, and 16.) The separation 
point, the separation shock angle, and the separation flow deflection angle were obtained 
from the schlieren studies. Separation and reattachment points were also measured by 
oil-flow studies. The trend in the movement of the separation point with a change in unit 
Reynolds number found from the oil-flow study clarifies the apparent contradictory results 
previously obtained over a much narrower range of Reynolds numbers (refs. 7 and 16). 
Reattachment pressures were determined from plots of the pressure distributions and the 
point of reattachment found from oil-flow studies. The model had a relatively short 
trailing-edge flap which undoubtedly had some effect on the extent of separation. 

The Lees and Reeves theory (ref. 24) was used to compute the pressure rise from 
the beginning of the interaction region to the end of the first pressure plateau region. The 
local similarity theory of reference 25 was used to calculate the upstream boundary-layer 
parameters used in the Lees and Reeves theory. The Lees and Reeves solution was 
joined to the upstream boundary-layer solution by matching the value of the momentum 
thickness at the beginning of the interaction region as determined from the experimental 
data. Detailed results of the calculations are presented herein. 

SYMBOLS 
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speed of sound; also velocity profile parameter, 


flow and 



for separated flow 
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specific heat at constant pressure 


specific heat at constant volume 
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stream function (see eqs. (6) and (7) and ref. 25) 


enthalpy; also used for heat-transfer coefficient in definition for Stanton 
number 


length of flat plate portion of model, 10 inches (0.254 meter) 


normalizing factor for theoretical solution, 1 foot (0.3048 meter) 


Mach number 


curve fit pressure gradient parameter 


Prandtl number 


Stanton number, h|cpp e u e 

6 t*/aU 

wall shear stress function, ==-(-= 
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static pressure 


unit free-stream Reynolds number per foot; also used in appendix B for 


dissipation function 
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total enthalpy function, 1 

h t,e 

temperature 

enthalpy function, E js w 5^* 

Stewartson's transformed velocity, (— )u 


velocity component parallel to surface 
coordinate parallel to surface 

Stewartson’s transformed coordinates (see ref. 24) 

velocity integral, ^ *^U^U e jdY 

pressure gradient parameter 

ratio of specific heats, Cp^c v 

trailing- edge flap angle 

boundary-layer thickness 

displacement thickness 


transformed boundary -layer thickness 

boundary-layer-displacement thickness ^ 11 (i - u/u e )ffi 


transformed displacement thickness, 6^* + e 
stagnation enthalpy ratio, h^h^g 
similarity variable (ref. 25) 


momentum thickness 



P 

P 


boundary- layer momentum thickness, ^ * S/E e (i - u/^X. 
mechanical "energy” thickness, ^ ’ l if ! 1 - H 2 /Ue% 


o u 

— e 


indicated shock angle (see sketch 1 on p. 14) 

viscosity 

density 

indicated flow deflection angle (see sketch 1 on p. 14) 


Subscripts : 

e local external conditions 

i transformed conditions 

o at beginning of interaction 

pi plateau value 

r reattachment value 

s separation point 

t stagnation conditions 

w wall conditions 

x along plate 

oo free- stream conditions behind oblique shock 


A prime indicates differentiation with respect to 77 . 


APPARATUS AND TEST PROCEDURES 


Description of the Models Used 

The pressure model used for the tests at the three wall temperature conditions is 
shown in figures 1 and 2. The model made from AISI type 347 stainless steel had a sharp 
leading edge about 0.001 in. thick (0.025 mm). The flat plate portion of the model was 

7.75 in. wide (0.197 m) and 10 in. long (0.254 m). A 2.0-in-long (0.0508 m) trailing- 
edge flap, which may be positioned at angles 6f of 0°, 10°, 20°, and 30° relative to the 
flat plate surface, extends across the back of the model. Tests were also made with 
upper side plates with sharp leading edges which extended back from the leading edge at 
an angle of approximately 6° as indicated by the line labeled "side plate" in figure 1. 

Instrumentation 

The pressure model was instrumented with 23 pressure orifices of 0.070-in. 

(1.78 mm) inside diameter as shown in figure 2. The instrumentation extends from 

4.75 in. (0.121 m) from the leading edge to 0.281 in. (0.00714 m) from the trailing edge 
of the flap. The pressures were measured with electrical hot-wire pressure gages and 
with electrical wire strain- gage type of pressure gages. The range of the electrical hot- 
wire gages is 0 to 20 mm Hg and the range of the strain- gage type varied from 0 to 1 to 0 
to 7.5 psia (0 to 6894.7 to 0 to 51 710 N/m2). The accuracy of the hot-wire pressure gage 
is thought to be ±0.05 mm Hg and the accuracy of the strain- gage pressure transducer is 
thought to be 0.75 percent of full-scale deflection. The calibration of the hot-wire gage 

is extremely nonlinear and was calibrated with a high degree of accuracy in the lower 
pressure range. The strain-gage transducers have a linear calibration over their rated 
pressure range. 

Pressure tests at T w /Tt = 0.43 and 0.14 were made with the electrical hot-wire 
pressure gages mounted inside the body of the model, as shown in figure 3. For pres- 
sure tests at T w ^Tt = 0.43 the strain-gage type was used also and was mounted inside 
the model injection box below the tunnel test section as shown in figure 4. Stainless- 
steel tubing and plastic vacuum tubing were used for all pressure leads. For the pres- 
sure tests at T w jTi = 0.74, both the electrical hot-wire and electrical wire strain-gage 
types were used with the gages mounted outside the tunnel, wrapped in a polyethylene bag, 
and immersed in an ice bath so that the gages could be kept at a constant temperature. 

Test Apparatus and Procedures 

The tests were conducted in the Langley Mach 8 variable -density hypersonic tunnel. 
The tunnel operated at a nominal Mach number of 8 over a range of Reynolds number per 
foot of 0.20 X 106 to 11.0 x 106 (per meter of 0.65 x 10® to 36.1 X 10*3). The tunnel 
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stagnation conditions varied from 25 psia and 1070° R (172.4 kN/m2 and 594° K) to 
2680 psia and 1470° R (18.5 MN/m2 and 817° K) for lowest and highest unit Reynolds 
numbers, respectively. A Mach number calibration of the facility can be found in refer- 
ence 26 and a further description is given in reference 27. Throughout the tests the 
model was set at 1/2° angle of attack resulting in a range of local Mach number on the 
plate from about 7.4 to 7.8. 

The pressure tests for T w /Tt = 0.43 were made with the model at essentially 
isothermal room temperature conditions since the data were taken at 1/2 sec after the 
model was positioned in the test section. Less than 1/2 sec was required for the pres- 
sures to reach an equilibrium value over the range of test conditions when the 0 to 20 mm 
gages were mounted inside the model. Equilibrium pressures were reached in approxi- 
mately the same time for the higher range pressure gages, which were mounted outside 
the model in the vacuum tight chamber because of their size. (See fig. 4.) 

The tests at T w /Tt = 0.14 were conducted with the model cooled by liquid nitrogen. 
The small hot-wire gages were mounted inside the model as shown in figure 3 and were 
sprayed with liquid nitrogen. The leading- edge piece had two passages drilled through 
it for liquid nitrogen cooling, as shown in figure 2. The pressure gage cavity (fig. 3) and 
back side of the flap were sprayed with liquid nitrogen. When the model reached an iso- 
thermal temperature of approximately 190° R (105.5° K) the model was injected into the 
tunnel and the pressure data were taken 1/2 sec after the model was positioned in the test 
section. The cold wall temperature of the model was measured by thermocouples mounted 
on the inside surface of the model at a point from the leading edge where the skin thick- 
ness was approximately 3/16 in. (0.00476 m) thick (fig. 2). A tight-fitting cover was 
placed over the surface of the model during the cool down period to prevent the formation 
of frost on the surface of the model. This cover was removed just prior to the injection 
of the model for a subsequent test. 

For the tests at T w /Tt = 0.74, the model was heated prior to exposure to the hyper- 
sonic stream to approximately 1000° R (555.5° K) by blowing subsonic heated air over 
the model. During an actual test when the model was in the hypersonic stream the model 
temperature settled out to a near equilibrium value which gave a wall-to-total tempera- 
ture ratio of approximately 0.74. The actual surface test temperature of the model 
ranged from 950° to 1200° R (528° to 666° K) for the lowest to the highest Reynolds num- 
bers, respectively. The thermocouples as shown in figure 2 were used to determine the 
surface temperature of the model for the hot wall test. The data were taken just prior 
to the time when there was a breakdown of the supersonic flow in the tunnel. The length 
of the run prior to flow breakdown was approximately 60 sec. 
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THEORETICAL APPROACH 


At the time the present tests were conducted, the most promising theoretical 
approach to the flat plate trailing- edge flap problem was the theory of Lees and Reeves 
in reference 24. (See appendix A.) Consequently, the method of Lees and Reeves was 
programed for use on the IBM 7090 electronic data processing system and approximately 
120 cases were calculated to check the agreement with experimental data. The first set 
of theoretical calculations was made for T w /Tt = 1.0, which was compared with the data 
for the adiabatic ^T w /Tt = 0.74) test conditions. Additional boundary -layer parameters 
were calculated for T w /Tt = 0.6 (appendix B) and compared with the data for room 
temperature wall conditions (T w /Tt = 0.43). 

The method of joining the Lees and Reeves separated flow solution to the upstream- 
boundary-layer solution as used herein deviates slightly from the method indicated in 
reference 24. The juncture of the two solutions is located at the beginning of the interac- 
tion region. (See fig. 5.) This point is determined from the results of experimental 
pressure distributions, rather than from an iterative solution from the complete theory. 
The beginning of the interaction region is defined as the point where the pressure first 
begins to rise above a flat plate value due to the adverse pressure gradient feeding 
upstream from the flap. In order to find the proper Lees and Reeves solution for a given 
experimental case the value of momentum thickness at start of the Lees and Reeves solu- 
tion was matched to the value of the momentum thickness from the upstream-boundary- 
layer solution at the x/L point defined as the beginning of the interaction. The 
upstream-boundary-layer calculation took into account the favorable pressure gradient 
due to the induced-boundary-layer effects. This method of matching was used because 
the complete Lees and Reeves solution from the beginning of the separation interaction 
region to the undisturbed flow downstream of the reattachment point could not be matched 
to the physical size of the model. Specifically, the calculations from the hinge line 
through reattachment and downstream to a point where the solution reached a Blasius 
type of flow yielded a streamwise x distance greater than the 2-in-long (0.0508 m) flap 
of the model. The calculations that were made with the Lees and Reeves theory were 
done with the intention of predicting the initial shape of the pressure curve and the level 
of the plateau pressure. For the initial 120 cases for T w /Tt = 1.0 and T w /Tt = 0.6 
only the region from the beginning of the interaction to the hinge line was calculated. 

The details of the Lees and Reeves calculations are given in appendixes A and B. 

Upstream-Boundary- Layer Calculation 

The local external flow properties for the upstream-boundary-layer solutions were 
calculated from the induced pressure effects for the case of weak interaction from the 
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theory of Bertram and Blackstock in reference 28. The equation for the induced pressure 
(for y — 1.4) is 


P _ 1 , a 

E, 01 1 

, 0.4285a 2 

7 , a 1 

\jl + a 

2(1 + a) 

(1 + a) 

2(1 + a)j 


where 

a = 0.7Gx 


G = 0.34416&L - 0.3859\ 

Vt ) 


— _ M 3 \/C 

A. 

C = 

M’oo'^'w 

The local values at the edge of the boundary layer of Mach number, temperature, velocity, 
and local Reynolds number per foot were computed from (y = 1.4). 


M e = ^P/Poo) 2 ^ ? ( 5 + M„ 2 | - 5.0 

(2) 

T e 5 + M^ 2 
T °° 5 + M e 2 

(3) 

u e = M e (49.1) ||T^ 

(4) 

1.812 X 10 8 (T e + 201.6)M e p e 
R e - 0 

rp Ct 

(5) 


The locally similar solutions of Beckwith and Cohen in reference 25 were used to 
calculate the upstream boundary layer with a pressure gradient as calculated by equa- 
tion (1). The boundary-layer equations in the similarity coordinates with the simplifying 
assumptions of constant Cp and pp and Np r = 1.0 reduce to 
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f ,M + ff M + p(% - f' 2 ) = o 


( 6 ) 


? n + fc = o 


( 7 ) 


where the notation is that of reference 25. 

Thus, with the local conditions external to the boundary layer calculated from equa- 
tions (1) to (5), the upstream boundary layer was obtained from a numerical solution to 
equations (6) and (7) at 35 points along the flat plate. 

Application of Theory to a Flat Plate 
With a Trailing- Edge Flap 

The Lees and Reeves solution between the beginning of the interaction region and 
the shock impingement point depends only on the separation point value of the transformed 
displacement thickness, the assumed local conditions at the edge of the boundary layer, 
and the previous upstream history of the boundary layer. Therefore, the Lees and Reeves 
solution, as calculated for a shock-wave— boundary- layer interaction, may be readily 
applied to a flat plate with a trailing- edge flap, from the beginning of the interaction region 
to the vicinity of the hinge line. 

The regions from the shock impingement point to beyond the reattachment point, for 
the case of a shock-wave — boundary -layer interaction, are assumed analogous to certain 
regions for a trailing-edge flap configuration if the wedge angle of the shock generator is 
half the flap angle since the final total compression angle of the inviscid flow would then 
be the same for the two configurations. The corresponding assumed flow models for the 
shock-wave — boundary-layer interaction and for the flat plate with trailing-edge flap are 
shown in figure 5 where it is assumed that the point of shock impingement and hinge line 
coincide. For the shock-wave — boundary-layer interaction a fluid element external to 
the boundary layer and moving parallel to the plate surface turns through an angle of 
8fj2 as it passes through the impinging shock wave. The same fluid element external to 
the boundary layer is turned again by an angle of 6f/2 as it passes through the reattach- 
ment compression fan and then moves downstream parallel to the plate surface. For the 
flat plate with a trailing-edge flap the fluid element external to the boundary layer is 
turned by an angle of 6f as it passes over the flap and through the reattachment com- 
pression fan. Thus, if the shock generator angle is taken as half the flap angle, the exter- 
nal flow for both conditions will experience the same total compression angle and will 
have a similar static-pressure history over the surface as is shown in part (c) of figure 5. 
In the flow models assumed for both configurations, the increase in entropy along the edge 
of the boundary layer is neglected - that is, the compression is assumed to occur through 
a series of weak waves rather than one strong shock wave. 
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TEST RESULTS 


Transition-Point Data 

It is known that transition can affect the extent of separation as well as the pressure 
levels associated with the region of separation; therefore, the location of the point of 
transition was determined from heat-transfer measurements made with a configuration 
which is geometrically similar to the pressure model of figures 1 and 2. (See ref. 29.) 

In figure 6 a typical example of the heat-transfer distribution obtained with the model 
tested as a flat plate (Sf = 0) illustrates where the position that was chosen as the point of 
transition is located (x/L « 0.90). 

The location of the point of transition for the 6f = 30° data in figure 6 is believed 
to be somewhere between the flat plate value of x/L ~ 0.90 and the point where the heat 
transfer in the separated region first begins to rise, x/L ~ 0.75. A summary of the flat 
plate transition data is shown in table I. 


TABLE I.- FLAT PLATE TRANSITION- POINT DATA 


Unit Reynolds number 

Distance from leading edge 

Local transition 

per foot 

per meter 

to beginning of transition, x/L 

Reynolds number 

2.65 X 106 

8.69 X 106 

0.90 

1.99 x 106 

3.45 

11.32 

.80 

2.30 

4.30 

14.11 

.75 

2.69 

10.90 

35.8 

.55 

5.00 


Over the range of unit Reynolds number an attempt was made to classify the type of 
separation as laminar, transition, or turbulent according to the definition of the three 
classes of separation as set forth in reference 20. The classification of the type of sepa- 
ration is readily defined at unit Reynolds numbers per foot of 2.65 x 10® (per meter of 
8.69 x 106 ) and above, where measurements of the flat plate transition point indicate a 
maximum downstream limit as to where transition should occur in the regions of separa- 
ted flow. For the separated flows at unit Reynolds numbers per foot below 2.65 x 10^ it 
is felt that transition could possibly occur in the separated region on the flap, somewhere 
before reattachment, even though there was no transition detected on the flat plate model. 
One possible reason for the occurrence of transition on the flap at lower unit Reynolds 
numbers is that the effect of separation and the effect of the flap deflection is to increase 
the local pressure and thereby increase the local Reynolds number in the separated 
region. Becker and Korycinski (ref. 21) found that in the presence of extensive separation 
the transition Reynolds number (based on free-stream conditions and distance to transi- 
tion point) was less by a factor of approximately 4 than the transition Reynolds number 
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for which there was no separation. The movement of the separation point with a change 
in unit Reynolds number (which will be discussed in detail in the section "Surface Oil- 
Flow Studies") also helps to classify the type of separation. 

Room Temperature Wall Pressure Tests 

Pressure distributions obtained on the plate with various flap angles at 
T w /T t = 0.43 for various Reynolds numbers are shown in figure 7. The measured pres- 
sures were divided by the pressures at the beginning of the interaction region. The pres- 
sures were obtained without side plates except for those presented in figure 7(d). All the 
pressure data in figure 7 are compared with the Lees and Reeves theory (ref. 24) for 
Tw/Tf = 0.6, and the beginning of the interaction region is determined as described in the 
section "Theoretical Approach." 

The agreement between experiment and the Lees and Reeves theory is, in general, 
good for the range of unit Reynolds number per foot from 0.22 X 10® to 1.46 x 10® (per 
meter from 0.72 X 10*3 to 4.79 x 106). The tests in this Reynolds number range are in a 
flow regime where plateau pressure level is not strongly affected by transitional effects. 
When transitional effects become more pronounced (Reynolds numbers per foot greater 
than 2.65 X 10® (per meter 8.69 x 10®)) the Lees and Reeves theory tends to underpredict 
the level of pressure in the area of the start of transition. The rise in pressure when 
transition occurs in the separated region was previously noted by Chapman, Kuehn and 
Larson in reference 20. For the tests at a unit Reynolds number per foot of 2.65 x 10® 
in figures 7(a) to 7(d) the rise in pressure above the laminar plateau value occurs at 
x/L « 0.90, which was previously given as the approximate transition location. For the 
tests at a unit Reynolds number per foot of 4.3 x 10® (per meter of 14.1 x 10^) in fig- 
ures 7(a) to 7(c), the start of transition is slightly upstream of the point of separation 
(based on the flat plate transition point); thus, little or no agreement would be expected 
at R = 4.3 X lO^/ft -with the laminar theory of Lees and Reeves. The side plate data in 
figure 7(d) at a unit Reynolds number per foot of 1.06 x 10® (per meter of 3.48 X 10^) 
show that as the flap angle increases the level of plateau pressure also increases, and, 
in addition, the Lees and Reeves theory tends to slightly underpredict for the 20° and 30° 
flap angles as was found for the model with no side plates. A comparison of the data for 
R = 2,65 x 10®/ft in figures 7(a) to (d) shows that, in general, the addition of side plates 
increases the extent of separation for all flap angles. The pressures measured upstream 
of the separation interaction region in figure 7 agree well with the viscous interaction 
theory of Bertram and Blackstock (ref. 28). The pressure ratio for the theoretical vis- 
cous interaction curves shown is divided by the theoretical pressure ratio at the beginning 
of the interaction. Flagged symbols in figure 7 indicate pressures measured off the center 
line of the model. (See fig. 2.) 
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In figure 8 are shown schlieren photographs for tests at T w /T t = 0.43 and flap 
angles of 10°, 20°, and 30° over a range of Reynolds number per foot from 0.22 x 10® to 
10.9 x 10® (per meter from 0.72 x 10® to 35.8 X 10®). These photographs were used to 
determine the angles of the leading-edge shock, the separated layer shock, and the deflec- 
tion angle of the separated boundary layer. In addition the photographs were used to 
determine the separation point based on the location where the separation shock wave 
intersects the boundary layer. These measurements of the separation point were found 
to be in good agreement with oil-flow separation- point data. 

Figure 9 shows the pressure distribution for a room temperature wall divided by 
the theoretical static pressure as calculated for the inviscid flat plate model set at 1/2° 
positive angle of attack for various flap angles and Reynolds numbers with and without 
side plates. The flap pressure ratios above values of 10 are not shown in figure 7 but 
are shown in figure 9. The separated-layer oblique-shock theory in figure 9 was calcu- 
lated from the apparent measured flow deflection angle above the separated layer (fig. 8) 
and based on the value of Mach number ahead of the interaction region. The value of 
Mach number ahead of the interaction region was calculated from the measured value of 
the leading-edge shock angle, above the interaction region, and the free-stream Mach num- 
ber. All oblique shock parameters were taken from reference 30. Also shown in fig- 
ure 9 are the values of the beginning of the interaction x 0 , the separation point x s , and 
the reattachment point x r . The values of the separation point and the reattachment point 
(fig. 9) were taken from oil-flow studies. The beginning of the interaction region is taken 
from the expanded plots, shown in figure 7, of pressures in the area of the plateau and 
interaction region, where the beginning of the interaction region is defined as previously 
noted in the section "Theoretical Approach." The separated-layer oblique-shock values 
of pressure agree, for the most part, with the measured pressure plateau values with the 
exception of a few of the lower unit Reynolds numbers where the shock pressures tend to 
slightly overpredict the measured values. 

The peak pressure ratio in figure 9(a) for the 30° flap at R = 10.9 x 106/ft 
(35.8 x 10®/m) is 29.6 compared with the inviscid-flap oblique -shock value of 28.6. The 
value of the pressure ratio reaches a maximum at approximately the same level as the 
inviscid-flap oblique-shock value apparently because of the thin boundary layer and the 
absence of any separation incurred at this highest Reynolds number. For 
R = 4.3 x 10®/ft and 2.65 x 106/ft (14.1 x 106/m and 8.69 x 106/ m ) the peak pressure 
ratios are approximately 42.0 and 45.0 as compared to inviscid values of 28.6 and 28.5, 
respectively. This peak value of about 45.0 for the lower Reynolds number represents 
an increase of 58 percent over the calculated inviscid value. The measured peak pres- 
sure exceeds the inviscid shock pressure because of the quasi-isentropic compression 
which occurs through a series of waves or several shocks rather than through a single 
shock. At Reynolds numbers per foot of 1.46 X 10®, 1.06 X 106, and 0.65 x 106 (per meter 
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of 4.79 x 106 , 3.48 x 106 , and 2.13 x 106), the pressure ratios have reached values greater 
than the corresponding inviscid values but apparently have not reached as high a value as 
they would have if the flap had been longer. This phenomenon is apparently due to the 
proximity of the reattachment point to the trailing edge of the flap, which does not allow 
the flow to turn completely parallel to the flap. This same effect becomes even more 
pronounced at the three lowest Reynolds numbers as the reattachment point moves closer 
to the trailing edge of the flap. In figure 9(b) a similar trend can be seen where for the 
20° flap, with the highest unit Reynolds number per foot of 10.9 x 106, the pressure 
approaches the inviscid value of 14.2, while at a Reynolds number per foot of 4.3 x 106 
the pressure rise increases to a peak and then starts an approach toward the inviscid 
value. In figure 9(c) for the 10° flap the peak pressure rise has a trend similar to that 
exhibited in figures 9(a) and 9(b) for the 30° and 20° flaps, only to a much lesser degree. 
The side plate data in figures 9(d) and 9(e), in general, indicate a slower rise to the peak 
pressure value on the flap. 

In figure 10 the indicated flow deflection angle xp (sketch 1), measured to within 
an accuracy of approximately ±0.3° from schlieren photographs, and the indicated shock 
angle 62 are plotted against unit Reynolds number. 



Sketch 1 

The angles were measured over the linear part of the shock wave and the separated 
region that was over the flat plate portion of the model. Figure 10 shows that for a given 
free-stream unit Reynolds number, the indicated flow deflection angle and corresponding 
indicated shock wave angle decrease as the flap angle decreases. Also, as the unit 
Reynolds number decreases the flow deflection and corresponding shock angle increase. 

A comparison of the plateau pressures obtained from oblique shock calculations, 
the Lees and Reeves theory, and measured static pressure values (for no side plates) is 
shown in figure 11. The oblique shock curves in figure 11 are faired curves taken from 
another plot of the individual data points. These faired oblique shock curves deviate 
approximately 2 percent from the individual data points and were calculated from the flow 
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deflection angles (fig. 10) by using the local values of Mach number just upstream of the 
beginning of the interaction. The local Mach number was calculated from the leading- edge 
shock- wave inclination and the free-stream Mach number. The Lees and Reeves predic- 
tion for the plateau pressure ratio for flap angles of 10°, 20°, and 30° was taken from 
plots of (p/p Q ^ j against M e?s for the various unit Reynolds numbers. (See appen- 
dix A.) The measured plateau pressures for flap angles of 20° and 30° were taken from 
the plotted data of figure 7. No experimental plateau pressures for a 10° flap angle are 
presented since the length of separated-flow region was so short that a pressure plateau 
could not form. In general, the oblique-shock value of plateau pressure is higher than 
the measured value and the Lees and Reeves value is slightly lower. (The difference in 
the Lees and Reeves value of plateau pressure between figure 7 and figure 11 lies in the 
fact that the plateau pressures for figure 11 are taken from plots of (p/P Q )pj against 
Me,s f° r the exact value of M e>s , whereas the theoretical plateau pressures in figure 7 
were taken from plots of p/p Q against x/L' for a nominal value of M e> s- (See appen- 
dix A.) 


Surface Oil- Flow Studies 

Technique for testing.- The pressure model (fig. 1) was used for the surface oil- 
flow studies and was tested with and without side plates. Prior to making an oil-flow test 
a pattern of drops consisting of an oil- lampblack mixture was placed on the surface of the 
plate as shown on the left in figure 12. The viscosity of the oil was increased with 
increasing unit Reynolds number in order that an oil-flow pattern could be established in 
approximately the same length of time for all runs. The model was rapidly injected into 
the wind-tunnel flow and when the apparent rearward movement of the oil drops stopped, 
giving a stabilized oil-flow pattern, the model was retracted from the stream. The tem- 
perature of the model in the area of separation varied from approximately 90° F to 100° F 
for all oil-flow tests . This temperature was measured with a thermocouple located 8 in. 
(0.2032 m) downstream from the leading edge on the inside surface of the model at an 
area where the skin thickness is approximately 3/16 in. (0.00476 m). The oil drops were 
placed on the model just prior to the run and a photograph was taken of the oil dots in the 
undisturbed position and then immediately after the run a second photograph was taken of 
the oil-flow pattern. These photographs were used as a means for determining the loca- 
tion of the separation point. A typical oil-flow pattern is shown in figure 12 for a free- 
stream unit Reynolds number per foot of 4.3 x 10® (per meter of 14.1 x 10®) and a flap 
angle of 20°. It can be seen from these oil -flow photographs that the surface streamlines 
are parallel to the inviscid flow before the interaction region; however, in the reverse 
flow region the surface oil flow shows the streamlines diverge indicating that the surface 
flow is three dimensional in nature. The use of the side plates considerably reduced the 
divergence of the flow in the separated region. 
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The viscosity of the oil for each Reynolds number and flap angle could not be chosen 
to insure a clear and distinct pattern of oil flow in the separated region for each run. For 
the majority of runs the oil droplets in the reverse flow region were displaced to some 
extent toward the separation point due to surface shear, where this movement could be 
detected by superimposing the negatives taken before and after the run. The use of the 
two negatives as an oil-flow evaluation technique was extremely helpful in determining 
the exact location of the point of separation which is indicated by the arrows. 

Separation poin t results .- The distance from the leading edge to the point of separa- 
tion as obtained from oil-flow data is shown in figure 13 for free-stream unit Reynolds 
numbers per foot from 0.22 x 10® to 4.3 X 106 (per meter from 0.72 x 10® to 14.1 x 106), 
for flap angles of 10°, 20°, and 30°, and for the studies with and without side plates. The 
results show that the separation point moves forward appreciably as the flap angle is 
increased and generally moves forward as the unit Reynolds number increases, up to 
R ~ 0.8 x 10®/ft (2.62 x 10®/m). At the lower Reynolds numbers (up to R = 0.4 x 106/ft 
(1.31 x 10®/m)) for the 10° and 20° flap angles, the rearward movement of the separation 
point with increasing Reynolds number is attributed to the fact that the outer edge of the 
viscous layer observed in the schlieren photographs passes well above the uppermost por- 
tion of the deflected flap (2 in. (0.0508 m) long). An examination of the schlieren photo- 
graphs showed that when the outer edge of the viscous layer was deflected in a direction 
parallel to the flap, the separation point began to move forward on the plate with increasing 
Reynolds number, as can be noted in figure 13 for the 30° deflected flap. It can be seen in 
figure 13 that for all three flap angles the separation point begins to move to the rear of 
the plate at a unit Reynolds number per foot of approximately 1.0 x 10® (per meter of 
3.28 x 106). This reversal in the trend of the separation- point movement is attributed to 
the flow going from a laminar separation to a transition type of separation. For these 
data it appears that separation up to a unit Reynolds number per foot of approximately 
1.0 x 10® gives laminar separation and for a unit Reynolds number per foot above approxi- 
mately 1.0 x 106 the separation is transitional, independent of the extent of separation and 
flap angle. The oil-flow tests were made with a model which had a trailing-edge flap 2 in. 
long. This short flap caused the reattachment point to be close to the trailing edge of the 
flap for the lowest unit Reynolds number tested. Undoubtedly the short length of the flap 
had some effect on the extent of the separation; however, to what degree the flap length 
affected the extent of separation would require a separate investigation. 

In figure 13 the effects of the side plates for all flap angles and Reynolds numbers 
is to move the separation point upstream with the exception of the 10° flap angle where for 
Reynolds numbers per foot from 0.65 x 10® (per meter from 2.13 x 106) to 4.3 x 10® 

(per meter from 14.10 x 106) the separation point moves slightly downstream with the 
addition of side plates. A comparison between the separation point, as indicated by the 
oil-flow pattern and as obtained from schlieren photographs, shows reasonably good 
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agreement when the schlieren separation point is chosen at the intersection of the shock 
wave from the separated layer and the apparent outer edge of the boundary layer. 

The effect of unit Reynolds number on the movement of the separation point.- From 
prior experimental investigations it was exhibited that when the Reynolds number was 
low enough the type of separation was purely laminar and that an increase in unit Reynolds 
number for any particular configuration moved the separation point upstream. For 
instance, the results from an experimental investigation by Miller, Hijman, and Childs at 
a Mach number of about 16 (ref. 15), in which all the data presented were classified as 
pure laminar separation, showed that as the Reynolds number increased (up to 
R ~ 1.0 x 10®/ft (3.28 x 106/m)) that the separation point moved upstream. The data of 
Chapman, Kuehn, and Larson (ref. 20) at Mach numbers of 2.7 to 3.5 show that, when the 
unit Reynolds number is increased so that the separation becomes transitional, the point 
of separation moves downstream with an increase in unit Reynolds number. In addition, 
the data show that a further increase in unit Reynolds number caused the separation to 
become turbulent and the separation point moved even further downstream. Becker and 
Korycinski (ref. 21) pointed out a similar effect on the movement of the point of separation 
on an axisymmetric flared body when the separation went from pure laminar to transi- 
tional and finally to a turbulent type of separation. The movement of the point of separa- 
tion with a change in unit Reynolds number shown in figure 13, and discussed in the pre- 
vious section "Separation point results," showed the same trend with a change in unit 
Reynolds number as found in references 15, 20, and 21. 

A set of theoretical calculations with the laminar theory of Lees and Reeves showed 
that an increase in the unit Reynolds number moved the separation point upstream. The 
calculation was made for a constant local Mach number at the point of separation (appen- 
dix A) and was based on the value of the displacement thickness found at the beginning of 
the Lees and Reeves interaction region and matched to the value of displacement thickness 
found on a flat plate with the local similarity solution of reference 25. 

Thus, the results from this paper and prior experimental investigations, as well as 
the theoretical calculations from the Lees and Reeves theory, show that for laminar sepa- 
ration the point of separation moves upstream with an increase in the unit Reynolds num- 
ber. Experimental results also show that for transitional and turbulent separation the 
separation point moves downstream with an increase in unit Reynolds number. 

Reattachmen t point results .- In figure 14 the results of an oil-flow study to deter- 
mine the reattachment point on the flap are shown. The reattachment point clearly moves 
toward the hinge line as the Reynolds number increases. The location of reattachment on 
the flap, for a given unit Reynolds number, varies only slightly with a change in the flap 
angle. It is only in the vicinity of a unit Reynolds number per foot of 1.06 x 10® (per 
meter of 3.48 x 106 ) that there is any appreciable variation in the point of reattachment 
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with a change in flap angle. This variation is attributed to a delay in the transitional 
effects with a decrease in flap angle. The solid lines in figure 14 are the predicted 
reattachment points as obtained from a linear extension of the dividing streamline from 
the point of separation (fig. 13) to the flap at an angle equal to the flow deflection angle of 
figure 10. The experimental data points are further downstream from the hinge line at 
the lower unit Reynolds numbers than the predicted values indicating that the dividing 
streamline curves upward as it moves downstream from the point of separation. 

In figure 15 the pressure ratios at the reattachment point on the flap are plotted 
against unit Reynolds number. The pressures in figure. 15 are limited in accuracy due 
to the difficulty in determining the precise point of reattachment in an area on the flap 
where the pressure gradient is often quite large. For a given unit Reynolds number the 
level of pressure increases with an increase in flap angle as might be expected; however, 
the reattachment pressure levels are considerably lower than the oblique shock inviscid 
pressure levels which are noted in figure 9. The reattachment pressure levels are much 
smaller than the inviscid oblique shock values indicating that the inviscid stream has not 
yet turned parallel to the flap and the compression process is still underway. 

Wall Temperature Effects 

In order to determine the wall temperature effect on separation, experimental tests 
were conducted at wall-to-total temperature ratios of 0.14, 0.43, and 0.74. Figure 16 
shows the measured flat plate static pressure, 5 in. (0.127 m) aft of the leading edge, 
divided by the theoretical inviscid flat plate static pressure at 1/2° positive angle of 
attack and the theoretical viscous interaction pressure calculated by the Bertram and 
Blackstock method (ref. 28) for various values of T w /Tt> The results in figure 16 show 
that the pressure ratios for T w /T^ = 0.43 and 0.74 decrease with an increase in unit 
Reynolds number and the theory overpredicts the measured pressure ratios by approxi- 
mately 10 percent. The pressure data at T w /Tt = 0.14 at the two lowest unit Reynolds 
numbers are considerably lower than the theory, whereas the remaining values of pres- 
sure ratio at the higher Reynolds numbers are lower than the theory by approximately 
10 percent or less. The large discrepancy between data and theory, for T w /Tt = 0.14, 
at the two lowest unit Reynolds numbers may be caused by a shift in the gage calibration 
due to a temperature effect on the hot-wire gage. 

Pressure paramet er d is tributi on for various flap angles.- The pressure tests con- 
ducted for various flap angles and wall-to-total temperature ratios are shown in figures 17 
to 22. The measured pressures have been divided by the measured pressure at the begin- 
ning of the interaction region. The pressure data in these figures for wall-to-total tem- 
perature ratios of 0.43 and 0.74 are compared to calculations based on the Lees and 
Reeves theory (ref. 24) for corresponding wall conditions of T w j Tj- = 0.6 and 
T w /Tt = 1.0, respectively. In addition pressure distributions based on the flat plate 
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viscous -inter action theory of Bertram and Blackstock (ref. 28) are calculated (and plotted) 
for wall-to-total temperature ratios of 0.43 and 0.74. The flagged symbols in figures 17 
to 19 indicate pressures measured off the center line of the model. (See fig. 2.) 

The data in figure 17 for 6f = 10° show reasonably good agreement with the theory 
of Lees and Reeves for unit Reynolds numbers per foot up to 1.46 X 10*3 (p er me ter up to 
4.79 X 10*3). The effect of wall cooling on the extent of separation and the level of plateau 
pressure is small for the three lowest Reynolds numbers. The data for the two highest 
Reynolds numbers show that, as the wall-to-total temperature ratio decreases, the extent 
of separation also decreases, particularly at the highest Reynolds number. 

An inspection of the experimental data of figure 17 shows that, as the wall-to-total 
temperature ratio decreases, the pressure gradient from the beginning of the interaction 
region to approximately the hinge line becomes slightly greater. Gadd, in reference 31, 
showed this effect of wall temperature on the pressure gradient rise for tests made at a 
Mach number of 3.0. Calculations with the theory of Lees and Reeves also show that, for 
T w /Tt = 0.6, larger pressure gradients occur in this region than occurred for 
T w /Tt = 1.0. In an earlier and different type of analysis from that of Lees and Reeves, 
Curie (ref. 32) and Gadd (ref. 33) showed that the pressure gradient at separation should 
be inversely proportional to the wall temperature, which is qualitatively the same finding 
that comes out of the Lees and Reeves solutions. The absolute value of the measured 
plateau pressure for T w /Tt = 0.74 is higher than for T w /Tt = 0.43. However, when 
the absolute values of plateau pressure are divided by the experimental value of the pres- 
sure at the beginning of the interaction, resulting values of the plateau pressure ratios for 
T w /T t = 0.43 and 0.74 are brought fairly close together. This trend can be seen from 
figures 16 to 22. 

The data of figures 20, 21, and 22 were obtained with side plates (fig. 1) at unit 
Reynolds numbers per foot of 1.06 X 10*3 and 2.65 X 10*3 (p er m eter of 3.48 x 106 and 
8.69 x 10*3). j n general, the results of figures 20, 21, and 22 show a slight decrease in 
the extent of separation with wall cooling at a unit Reynolds number per foot of 1.06 x 10*3 
and a marked decrease in the extent of separation with wall cooling at a unit Reynolds 
number per foot of 2.65 x 10*3. Figures 20 to 22 also show a decrease in flap pressure 
ratio with an increase in the wall-to-total temperature ratio for both unit Reynolds 
numbers. 

Separation point data .- A summary of the separation point data taken from the 
schlieren photographs of figure 23 for various wall-to-total temperature ratios and two 
flap angles are shown in figure 24. The schlieren photographs were taken during the 
pressure tests shown in figures 18 and 19. For the 30° flap angle, the results in figure 24 
for the room temperature wall (T w jT^ = 0.43) and the hot wall (T w /Tj. = 0.74) show that 
the separation point moves forward on the plate with increasing unit Reynolds number per 
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foot up to a value of approximately 1.0 X 10$ (per meter, 3.28 X 10®). An increase in 
unit Reynolds number per foot above approximately 1.0 x 10® for all three wall-to-total 
temperature ratios rapidly moves the separation point toward the hinge line for both the 
20° and 30° flaps. For pure laminar separation the wall temperature effects for a flap 
angle of 30° show that the highly cooled wall (T w /Tt = 0.14) has the greatest extent of 
separation, whereas for T w /Tt = 0.43 and 0.74 nearly an equal extent of separation is 
noted. The wall temperature effect for 6f = 20° is quite similar to that for the 30° flap 
with the extent of laminar separation being the greatest for the highly cooled wall 
(T w /Tt = 0.15), the room temperature wall T^ = 0.43) having the smallest extent of 
separation, and the hot wall = 0.74) having an extent of separation falling between 

the other two wall conditions. For both the 20° and 30° flap angles, after transition 
occurs the effect of wall cooling for a given unit Reynolds number tends to decrease the 
extent of separation significantly. (Transition is assumed to occur near a unit Reynolds 
number per foot of approximately 1.0 X 10® based on previously discussed results at 
T w /Tt = 0.43.) This wall temperature effect, decreasing the extent of separation, is par- 
ticularly noticeable at unit Reynolds numbers per foot above 2.65 x 106 (per meter above 
8.69 x 10®). The effect of wall cooling in reducing the extent of separation was also 
shown in reference 34 for a cone- cylinder flare at Mach number 5 over a similar unit 
Reynolds number range. 

Plateau pressur e data.- In figure 25 the values of plateau pressure ratios for various 
wall-to-total temperature ratios and flap angles were taken from the plots of figures 18 
and 19. The level of plateau pressure decreases with an increase in the unit Reynolds 
number for all three temperature ratios and both flap angles. The highly cooled wall 
(T w /Tt = 0.14) and the hot wall (T w j T^ = 0.74) give higher values of plateau pressure than 
the room temperature wall (T w /Tt = 0.43) up to a unit Reynolds number per foot of 
1.0 x 10® (per meter of 3.28 x 10®). At a unit Reynolds number per foot of about 1.0 X 106 
the transitional effects become significant and the plateau pressure level, for a given unit 
Reynolds number, decreases with an increase in wall cooling. 

Flow deflection and shock wave angles.- The schlieren photographs (fig. 23) were 
used to determine the wall temperature effect upon the flow deflection angles and the shock 
wave angles shown in figure 26. An analysis of figure 26 shows that the flow deflection 
angle and the shock wave angle, measured at the point of separation, decrease continuously 
with increasing unit Reynolds number for all three wall-to-total temperature ratios. The 
effect of wall temperature on flow deflection and shock angle is not large; however, it 
should be noted that the highly cooled wall (T w /Tt = 0.14) and the hot wall (T w ) Tt = 0.74) 
both have slightly higher values of deflection and shock wave angle than does the room 
temperature wall ^T w / T^ = 0.43). This fact tends to confirm the observations noted pre- 
viously regarding effect of cooling on the measured pressures. 
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Plateau pressure correlation .- Figure 27 shows a correlation of the plateau pres- 
sure coefficients in terms of the local Reynolds number and local Mach number at the 
beginning of the interaction region. The local flow properties were evaluated from the 
measured wall pressures, shock angles, and the oblique shock relations. These data 
show reasonable agreement with an expression given by Hakkinen, Greber, Trilling, and 
Abarbanel in reference 35 but fall slightly above the theory of Erdos and Pallone in ref- 
erence 36. The lower Mach number data of Chapman, Kuehn, and Larson in reference 20 

and the data of Miller, Hijman, and Childs in reference 15 also are higher than the expres- 

1/4 

sion for (cp)p](Re,x) 0 °f Erdos and Pallone. 

In figure 28 the plateau pressure data of the present investigation, as shown in fig- 
ure 27, are compared with the data of references 4, 11, 15, 19, 20, 23, and 35. The corre- 
lation of (cp)pj(Re,x)o"^ is plotted as a function of the local Mach number at the begin- 
ning of the interaction Mg }0 from 1.3 to 14.3. Even though there is scatter in the data 
the agreement with the theories of references 35 and 36 can be seen over the range of 
Mg 0 . In particular the agreement with the theory of reference 35 appears reasonably 
good up to a Mach number of approximately 8, whereas for Mg )0 values greater than 8 
the data show closer agreement with the theory of reference 36. The curves obtained 
in references 35 and 36 are somewhat dependent upon the assumed relationship between 
Reynolds number and skin friction. The use of other methods of determining skin fric- 
tion could lead to other curves than those shown. 

CONCLUDING REMARKS 

An experimental investigation of pressures on a flat plate with various flap angles 
and three wall-to-total temperature ratios was made. The model tested had a short 
trailing edge flap (2 in. (0.127 m) long) which, for the lower unit Reynolds numbers, 
resulted in the point of reattachment being near the trailing edge of the flap. The short 
length of the flap and the associated reattachment point close to the trailing edge of the 
flap undoubtedly had some effect on the extent of separation and clearly had some effect 
on the maximum pressure rise on the flap. 

Experimental pressures in the separation region were compared with calculations 
made by the Lees and Reeves theory (AIAA Journal, Nov. 1964) for a shock-wave — 
boundary- layer interaction. Only the calculation for the region from the beginning of the 
interaction region to the maximum extent of the laminar plateau region was used due to 
the short length of the flap. Agreement between the separation pressure data and the 
Lees and Reeves theory was found to be reasonably good. Lees and Reeves upper and 
lower branch boundary -layer parameters for a wall-to-total temperature ratio of 0.6 
were calculated from local similarity boundary -layer solutions for use in the separated 
flow solution. 
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The results of an extensive surface oil-flow study conducted at room temperature 
wall conditions showed that for pure laminar separation an increase in unit Reynolds num- 
ber moved the separation point upstream; however, when the separation becomes transi- 
tional, an increase in unit Reynolds number moved the separation point toward the hinge 
line. It was also found that the use of side plates and an increase in the size of the flap 
angle increased the extent of separation. The variation of the extent of separation with a 
change in unit Reynolds number was nearly the same for conditions with and without side 
plates. The flap reattachment point, according to the results of the oil-flow study, showed 
that for a given unit Reynolds number the point of reattachment changed very little with 
a change in flap angle from 10° to 30°. However, the pressures at reattachment varied 
considerably. 

The maximum pressures measured on any particular flap varied considerably with 
a change in unit Reynolds number. The maximum pressures measured on the 30° flap at 
a unit Reynolds number per foot of 2.65 X 10^ (per meter of 8.69 X 106) indicated a pres- 
sure rise of approximately 58 percent over the oblique shock reattachment value. The 
reason for this rise in pressure is believed to be due to a quasi-isentropic compression 
occurring through a series of waves (or shocks) rather than a single shock. Similar 
pressure rises, only to a lesser degree, are noted for the 20° flap angles for similar unit 
Reynolds numbers. 

The measured separation flow deflection angles and shock-wave angles for wall-to- 
total temperature ratios of 0.14, 0.43, and 0.74 indicate qualitative agreement with calcu- 
lations made from measured pressures and the theoretical calculations made from the 
Lees and Reeves theory. 

A change in the wall-to-total temperature ratios (T w /Tt) showed: 

(1) For the range of unit Reynolds number that appeared to give laminar separation, 
T w /Tt ~ 0.14 gave the highest plateau pressure and the greatest extent of separation; the 
smallest value of plateau pressure and the least extent of separation occurred for 
Tw/Tf ~ 0.43; and for T w /Tt ~ 0.74 the value of plateau pressure and length of separa- 
tion was between these two extremes. 

(2) For the range of Reynolds numbers where the separation was classified as first 
becoming transitional (with transition believed to occur in the separated layer very near 
the flap) the level of plateau pressure and the extent of separation were in reasonably 
close agreement for all three flap angles. 

(3) For the highest Reynolds number, where the separation was clearly transitional 
(with transition occurring in the separated layer above the flat plate portion of the model), 
a decrease in the wall-to-total temperature ratio markedly decreased the extent of sepa- 
ration and also reduced the level of the plateau pressure. 
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The viscous interaction theory of Bertram and Blackstock (NASA Technical 
Note D-798) agreed well with the experimental flat plate data upstream of the separation 
effects for T w /Tj. = 0.43 and 0.74. 

A correlation of the plateau pressure with the local Mach number and Reynolds num- 
ber at the beginning of the interaction showed reasonable agreement with the theory of 
Hakkinen, Greber, Trilling, and Abarbanel (NASA Memorandum 2-18-59W). 

The effect of side plates on the model as compared to the results with no side plates 
showed a reduction in the three dimensionality of the flow in the separated region, a slight 
increase in the plateau pressure level, and a considerable increase in the extent of 
separation. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 29, 1967, 

129-01-08-38-23. 
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APPENDIX A 


SOLUTION FOR THE LAMINAR SEPARATED BOUNDARY LAYER 
Review of Theoretical Literature 

The problem of flow separation has been investigated theoretically since the time of 
Prandtl’s early works published in reference 37. In recent times one of the first efforts 
toward an analysis of supersonic separation was made by Chapman, Kuehn, and Larson 
(ref. 20) in which the results of Chapman's mixing layer analysis (ref. 38) were used. 
Chapman’s analysis represents a limiting case for separation with the assumption that 
the boundary -layer thickness is zero at the point of separation. 

The Karman-Pohlhausen method was used by Gadd, Curie, and Savage (refs. 39, 32, 
and 40, respectively) without a great deal of success primarily because the assumed 
velocity profiles in the region of separation did not give the reverse flow found in experi- 
ment. Crocco and Lees, in reference 41, developed a semiempirical method which 
depends on the rate of entrainment of fluid from the external stream into the boundary 
layer. Results from the Crocco- Lees method were only in qualitative agreement with 
experimental data. The Crocco-Lees method, modified by Glick (ref. 42), predicted 
results that were in good agreement with pressure data as obtained from experiment. 

This method uses the concept of the dividing streamline; however, empirical data are 
required for its application. The Crocco-Lees method was also used by Bray, Gadd, and 
Woodger (ref. 43) with reasonable success. 

Tani in reference 44 used an analysis similar to that of Wieghardt (ref. 45) and 
Walz (ref. 46) in that his solution for an attached flow with an adverse pressure gradient 
used the first moment of momentum in addition to the zeroth moment and continuity equa- 
tions. Tani used a quartic representation for the velocity profiles; however, the boundary 
condition which required that the momentum equation at the wall be satisfied was dropped. 
When this boundary condition was neglected the resulting one parameter, which charac- 
terizes the family of velocity profiles, was not directly related to the static pressure dis- 
tribution but was directly proportional to the shear stress at the wall. The one parameter 
describing the family of velocity profiles was obtained from the simultaneous solution of 
the zeroth moment of momentum and the first moment of momentum equations. The 
results of Tani's analysis have been found to be in good agreement with exact solutions of 
the boundary -layer equations. Poots, in reference 47, extended Tani’s method by adding 
the energy equation to the continuity and two momentum equations. 

Abbot, Holt, and Nielsen, in reference 48, studied the separated flow problem by 
using the continuity equation, the zeroth and the first moment of momentum equation, and 
the energy equation with a fourth degree polynomial expression for the velocity and 
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temperature profiles and with one undetermined parameter per profile. The resulting 
separated flow pressure distributions did not have the correct trends primarily because 
of the use of polynomials for the velocity and temperature profiles. Lees and Reeves, in 
reference 24, developed a method for the shock-wave — boundary-layer interaction prob- 
lem wherein the continuity equation, the momentum equation, and the first moment of 
momentum equations are solved simultaneously with a one parameter family of velocity 
and enthalpy profiles. Lees and Reeves used the Cohen and Reshotko profiles for the 
highly cooled wall cases and added the Stewartson profiles for the adiabatic wall, as found 
in references 49 and 50, respectively. This method gives good agreement with experi- 
mental pressure data obtained in the present paper for both adiabatic and cooled walls 
I T^- = 0.6); however, for the highly cooled wall (T w jT^ = 0.2) and for quantitative heat- 
transfer predictions, indications are that the method is inadequate. The most promising 
method for predicting both pressure and heat transfer under highly cooled wall conditions 
seems to be that of Holden (ref. 51), which adds the energy equation to the conservation of 
mass and the zeroth and first moment of momentum equations. Holden's method of solu- 
tion is similar to that of Lees and Reeves in that he uses the velocity and enthalpy profiles 
from the upper and lower branches of the Cohen and Reshotko solution (ref. 49). However, 
Holden's family of velocity and enthalpy profiles are determined by two parameters, one 
of which defines the velocity profile and the other defines the enthalpy profile. Both 
Holden and Lees and Reeves in their methods of solution uncouple the boundary-layer 
velocity profiles from the pressure gradient parameter associated with the Cohen and 
Reshotko solution. For the Lees and Reeves method of solution once the velocity profile 
is determined there is only one enthalpy profile associated with the given velocity profile. 
On the other hand in Holden's method the enthalpy profile is uncoupled from both the pres- 
sure gradient parameter and the velocity profile, and with the inclusion of the energy equa- 
tion the enthalpy profile parameter can be determined. The results of Holden's method 
agree well with his highly cooled wall experimental heat-transfer and pressure data 
(ref. 51). 


Lees and Reeves General Method of Solution 

The Lees and Reeves theory (ref. 24) gives the solution for the laminar boundary 
layer in which a pressure disturbance is propagated upstream through a supersonic flow. 

A pressure disturbance may be generated by a shock wave impinging on the boundary layer 
or it may be caused by a trailing- edge flap - the Lees and Reeves theory may be applied 
to either of the two types of disturbances. (See fig. 5.) The theory, and its calculations 
as used in this paper, applies at the beginning of the interaction region and is used to the 
point of the shock impingement (the hinge line) for T w /Tt = 0.6 and 1.0. The flow region 
on the flap was calculated for a few cases but it was found that the extent of the calculated 
region over the flap was considerably longer than the actual size of the flap on the model. 
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The Lees and Reeves method requires the solution of the previously mentioned con- 
servation equations coupled with an inviscid-streamline — Prandtl-Meyer solution. The 
method gives the solution for the boundary layer and the flow external to the boundary 
layer within the framework of a single parameter family of velocity profiles. This one 
parameter determines the velocity and enthalpy profile for specified regions of attached 
and separated flows and yet is not directly related to the local static-pressure gradient. 
The local pressure gradient is determined from the local inviscid flow inclination and the 
Prandtl-Meyer solution. 

As previously noted for T w /Tt = 1.0, the Stewartson (ref. 50) and the Cohen and 
Reshotko (ref. 49) boundary- layer solutions were used to evaluate the integral parameters 
used in the Lees and Reeves solution. With heat transfer, T w /Tt = 0.6, the boundary- 
layer profiles were calculated from the local similarity solutions with Np r = 1.0, 

Cp = Constant, and p\i = Constant. These boundary- layer solutions were used to calcu- 
late the integral parameters. (See appendix B.) These integral parameters were then 
curve fitted (by a polynomial expression) as a function of the single parameter a used 
to describe the entire family of velocity profiles for both attached and separated flow. 

The coefficients of the polynomials from the fitted curves of the integral parameters for 
Tw/Tf = 0.6 are tabulated in tables II and in. The local similarity solutions reduce to 
a solution which is the same as that of Cohen and Reshotko (ref. 49) when a Prandtl num- 
ber of unity and a constant heat capacity are used. 

The Lees and Reeves method of calculation starts at the point of separation and 
moves upstream until a flat plate (Blasius type) solution is reached at the upstream end of 
the interaction. The values of local Mach number and unit Reynolds number are fixed at 
the point of separation and the value of the transformed displacement thickness at the 
point of separation is iterated for until the proper upstream solution is found at the begin- 
ning of the interaction region. The conditions needed for satisfying the two-point boundary 
value problem are the rate of change with respect to the transformed X distance of the 
local Mach number and the shape parameter approaching zero at the same time as the 
parameter describing the family of velocity profiles approaches the zero pressure gradi- 
ent value for attached flow. After the correct value of the displacement thickness, at the 
point of separation, is found the solution moves downstream into the separated flow region. 
The basic equations used in the Lees and Reeves method of solution are found in refer- 
ence 24 and were integrated by a fourth-order Runge-Kutta integration procedure which 
extrapolates to a zero interval size as a correction factor. The equations were all inte- 
grated with 6t* as the independent variable. 
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Procedure Used for Application of the Lees and Reeves Theory 

In figure 29 the experimental values are shown of the beginning of the interaction 
region as taken from measured pressure distributions at three different wall-to-total tem- 
perature ratios. The beginning of the interaction region (x/L) Q is selected as the point 
where the pressure begins to rise above the undisturbed upstream values due to the 
adverse pressure gradient feeding forward from the flap. The Lees and Reeves solution 
is then joined to the upstream boundary- layer solution at this (x/L) Q . The upstream 
boundary- layer momentum thickness is shown in figures 30 and 31 for a plate 10 in. 

(0.254 m) long at T w /T{ = 1.0 and 0.6, respectively. The calculations in figures 30 and 
31 were made for a unit Reynolds number range per foot from 0.22 x 10® to 4.3 x 10® (per 
meter from 0.72 x 10® to 14.11 x 10®), respectively. (The Mach number and unit Reynolds 
number actually varied slightly along the plate according to the weak interaction equa- 
tions.) The momentum thickness at the beginning of the interaction from the Lees 

and Reeves solution for the various test Reynolds numbers is plotted against the assumed 
local Mach number at the point of separation M e s in figures 32 and 33. Specific cases 
using the Lees and Reeves theory were calculated for unit Reynolds numbers per foot 
from 0.22 x 106 to 4.3 x 10*5 and local Mach numbers at separation from 6.5 to 7.5 at 
0.1 intervals in M e g . The undisturbed flat plate Mach number was varied from 7.4 to 
7.8 to correspond to the change in the test section Mach number with a change in free- 
stream unit Reynolds number (ref. 27). Typical plots of pressure against x/L' from 
the beginning of the interaction region are shown in figure 34 for various values of unit 
Reynolds number and T W ^T^ = 1.0 and 0.6. 

The first step in the application of the Lees and Reeves theory as used in the pres- 
ent paper is to obtain from figure 29, for a given T w /T(. and unit Reynolds number, the 
(x/L) q value for the beginning of the interaction. This value of (x/L) Q is then used to 
find the value of the momentum thickness at the beginning of the interaction 0 Q for the 
upstream solutions from figure 30 or 31. This value of 0 o is then used in either fig- 
ure 32 or 33 for a given value of T w /Tt and R to obtain the corresponding value of 
M e>s , which in turn is used to specify the particular Lees and Reeves solution and pres- 
sure distribution (fig. 34) for the given test conditions. 

It can be seen in figure 34, for all values of unit Reynolds number, that as the Mach 
number at separation increases, the rate of the pressure increase (with distance) 
decreases for a sizable x/L* distance before the pressure begins to climb toward a 
plateau value. This slow rate of increase in pressure occurs in the region between 
pjp Q = 1-00 and 1.10 and makes it difficult to determine the beginning of the interaction 
region. For the purpose of comparing the theory with experimental data this difficulty 
was overcome by linearly extrapolating the slope of the curve at P/P Q =1.10 to the 
abscissa as is shown typically by the dash-dot line in figure 34(a) for R = 0.22 x 10®/ft 
(0.72 x 106/m) and M e s = 6.8. The point where the extrapolated lines crossed the 
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x/L' abscissa was considered to be the theoretical point of the beginning of the interac- 
tion region; this was matched to the experimental value of the beginning of the interaction 
region. Thus, the theoretical and experimental techniques of determining the beginning 
of the interaction were consistent, in that both methods used a sudden pressure rise to 
define the beginning of the interaction region. 

It should be noted in figure 34 that, at the separation point for T w /Tt = 0.6, the 
pressure-ratio curves have a distinct discontinuity in slope. When the theoretical curves 
were compared with the experimental data the discontinuities in slope were faired to give 
a smooth pressure rise. The reason for the discontinuity in the slope of the pressure 
curves at the point of separation is believed to be due to the change in slope of the sepa- 
rated and attached profile parameters upstream and downstream of the point of separation. 
(See appendix B.) For the majority of the cases calculated in figure 34 the plateau pres- 
sure level for the adiabatic wall conditions (T w )Tt = 1.0) was slightly higher than that for 
the cool wall (T w j T^. = 0.6), where the extent of separation was long enough to allow the 
pressure ratios for both wall conditions to reach an almost constant plateau value. 

In figure 35 a comparison is shown between the growth of the momentum and dis- 
placement thicknesses for the upstream similar solution boundary layer (T^jT^ = 0.43) 
and for the downstream Lees and Reeves solution at a unit Reynolds number per foot of 
0.22 x 106 (per meter of 0.72 x 106 ), where the two solutions were joined based on 
(x/L) Q for 6f = 30° and M e s = 6.6. The results of a typical calculation shown in fig- 
ure 35 indicate that the momentum thickness changes less than the displacement thickness 
in the presence of an adverse pressure gradient for the interaction and separated regions. 

Plots of the plateau pressure against the local Mach number at separation as 
obtained from the Lees and Reeves theory are shown for T w /Tj- = 0.6 and 1.0 in fig- 
ures 36 and 37, respectively. The dashed curves represent the theoretical values of 
plateau pressure based on the experimental value of the beginning of the interaction found 
in figure 29. In addition to figure 29, figures 30 to 33 were used to obtain the proper 
value of the local Mach number at separation (as discussed previously) for the dashed 
curves in figures 36 and 37. The predicted values of plateau pressures in figures 35 and 
36 are confined to a narrow band which decreases steadily with an increase in unit 
Reynolds number. 


28 



APPENDIX B 


LEES AND REEVES BOUNDARY- LAYER PARAMETERS FOR T w /Tt = 0.6 


The curve fit parameters for the boundary- layer integral parameters for 
T w /T t = 0.6 that were used in the Lees and Reeves calculation are given in tables II and 
HI for attached and separated flow parameters, respectively. These parameters were 
calculated from local similarity solutions using a Prandtl number of unity and constant 
heat capacity. The integral parameters, in the notation of reference 25, are 


p6 

f’(l - f’)d?7 

.^0 

p6 

(S-f’)d77 
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C 6 ft(i . f' 2 )drj 
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The value of f’, at the edge of the boundary layer, was taken as 0.9995 for all calculations. 
The parameters in equations (Bl) to (B6) were fitted to various order polynomials by the 
method of least squares as a function of the parameter a defined by the Lees and 
Reeves method. The parameters H, J, P, R, T*, and Z in tables n and in are 
defined in equations (Bl) to (B6) while and m^ are defined in reference 29 where 
the expression for heat transfer is given. 

The polynomials in tables n and in occur in the general form 

H = A + Ba + Ca2 + Da3 + Ea 4 + Fa 5 + Ga6 + Ha? 

and vary from fourth order for the attached flow parameter to as high as seventh order 
for separated flow. 


TABLE H.- COEFFICIENTS FOR ATTACHED FLOW 

T w /T t = 0.6 


Parameter 

A 

B 

C 

D 

E 

H 

0.3321009 

0.13986596 

0.01179970 

-0.00946576 

0.00275081 

J 

.50348534 

.19925458 

.02442816 

.00089497 

-.00135237 

R 

.9388926 

-.41635211 

.19676828 

-.06663378 

.01160929 

X P (if a > 0.4) 

.00432494 

.28136774 

-.04773970 

-.01147897 

.00336312 

Z 

1.894155 

.79107769 

.18536477 

-.00701481 

.007183026 

T * 

.17239089 

.13440415 

.01340911 

.01577639 

-.00074312 

?w 

.12409628 

.05559569 

-.03955804 

.02479668 

-.00555449 

m „3 

-.1094301 

-.01730765 

.05114864 

.02927415 

.01006094 


% 0 Sag 0.4, P = 0.2715a. 


TABLE m.- COEFFICIENTS FOR SEPARATED FLOW 
T w /T t - 0.6 


Parameter 

A 

B 

C 

D 

E 

F 

G 

H 


H 

0.3318532 

-0.52778371 

0.2626242 

-3.0420035 

4.3360826 





J 

.50378811 

-.99282899 

5.100377 

-52.405373 

251.93463 

-630.63668 

788.67751 

-384.52197 

*R (if 

a > 0.06) 

.92651705 

3.5312703 

-48.027304 

502.05092 

-2393.2702 

6017.3521 

-7518.2889 

3648.4459 

2 P (if 

a > 0.05) 

-.00302374 

-.5727471 

-7.2003253 

53.081567 

-263.80331 

669.6665 

-745.74645 

286.17898 


Z 

1.8961545 

! -2.5925828 

-6.8986577 

55.663569 

-186.92544 

281.81194 

-155.92 



T* 

.17241606 

-.40731561 

-.11445041 

.8055802 

-1.1635996 

1.1703974 



3 C {ji 

a > 0.05) 

.12554626 

-.16814151 

1.4798979 

-19.391634 

105.05788 

-293.36703 

401.23095 

-210.28584 


m /3 

-.10887434 

-.10113126 

2.5761168 

-17.154012 

59.001127 

-89.996992 

48.966367 



ilf (0 £ a ^ 0.06), R = 0.9388926 + 1.806118a. 

2 If (0 ^ a £ 0.05), P = -0.889a. 

%(0^aS 0.05), « 0.12409628 + 0.1022052a. 
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Figure 5.- Comparison of the basic flow models assumed for a shock-wave— boundary- layer interaction and flat plate with a 
trailing-edge flap and their associated laminar pressure distribution. 
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The effect of transition on heat transfer for a flat plate = 0°) and for a deflected flap ( = 30°). 
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(q) R - 10.9 x 10% (35.8 x 10%). 


(z) R - 10.9 x 10% (35. 8 x 10%). 


L-67-6601 

Figure 8.- Schlieren photographs of the flow separation model at T w ^ = 0.43 for three flap angles and various unit Reynolds numbers. 
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Figure 18.- Concluded. 







Figure 19.- Continued, 
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Figure 23.- Concluded. 
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Figure 26.- The effect of Reynolds number, wall temperature, and flap angle on the indicated separation point flow deflection angle and the 

indicated separation shock wave angle. 
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Figure 28.- A correlation of plateau pressure with Reynolds number and Mach number. 
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Figure 34.- The theoretical pressure rise due to separation for various local Mach numbers at separation for cool and adiabatic wall conditions at various 

unit Reynolds numbers. 
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Figure 36.- The effect of Reynolds number and flap angle on the theoretical prediction of plateau pressure for cool walls. 
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